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Abstract 

Using a d x 2_ y 2 superconductor in 2+1 dimensions we show that the Nambu Goldstone fluctu- 
ations are replaced by dissipative excitations. We find that the nodal quasi-particles damping is 
caused by the strong dissipative excitations near the nodal points. As a result we find that the 
scattering rates are linear in frequency and not cubic as predicted in the literature for the "d" 
wave superconductors. Our results explain the recent angle resolved photoemission spectroscopy 
and optical conductivity in the BSCCO high T c compounds. 

PACS numbers: 
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The recent angle-resolved photoemission spectroscopy (ARPES) [1] and optical conduc- 
tivity [2] studies of the superconducting high T c compound BSCCO show scattering rates 
which are linear in temperature [1,2,3]. Moving away from the nodal direction the scattering 
rates appear to level of and become temperature independent. These scattering rates are in 
strong contrast to the uj 3 scattering suggested for the d-wave like pairing [4]. 

Our results are expected to hold for any d x 2_ y 2 superconductor in 2+1 dimension. The fact 
that our predictions are only observed in BSCCO suggest that the YBCO superconductors 
have additional d-wave order parameters u id xy \ "s" or have a large hopping in the "c" 
direction. In both situations the dissipative behavior caused by massless Dirac fermions in 
2+1 is absent! 

Quasi-particles properties of nodal fermions in d x 2_ y 2 superconductors and 2+1 dimen- 
sions are investigated. We find that the Nambu-Goldstone fluctuations or the massive plasma 
mode caused by the Coulomb interaction is replaced by a dissipative collective mode - a gauge 
invariant field. This result follows from the fact that the Nodal-Dirac [5-9] fermions which 
couples to a gauge field in 2+1 dimension induces critical Q.E.D. and not regular photons 
[5-10]. We find that this dissipative mode causes the quasi-particle self energy in the vicinity 
of the nodal points to have scattering rate which are linear in frequency [1-3] . Moving away 
from the nodal point one finds that the scattering rate is reduced. When the d-wave order 
parameter has additional components, u id xy " , "s", a' or single electrons hopping occurs in 
the "c" direction the scattering rate of the quasi-particles is suppressed. 

In order to understand the origin of the dissipative mode we follow S. Weinberg [10] and 
replace the fermion operator C a (x) by a neutral fermion field C a (x) = exp(—ia(x)/2)C a (x) 
where "a" is the Nambu-Goldstone phase. In the presence of an electromagnetic field a ext 
one finds [10] that the electromagnetic response is represented in terms of the gauge invariant 
fields a M = {d^a — 2a^*). For any superconductor the electromagnetic response contains two 
parts, the diamagnetic contribution given by the vector potential square, \p{da — 2a ext ) 2 (p 
is the electron density) and the paramagnetic response which is obtained by second order 
perturbation theory. 

For an "s" wave superconductor the density-density response function is proportional 
to the superconducting density [12]. The current-current response function represents the 
parametric polarization, T^'^ which at T = and q — > vanishes, = 0! Therefore the 
diamagnetic term \p{da — 2a ext ) 2 is not normalized. As a result the Lagrangian takes the 
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form; L sc ~ \p ^i(d t a — 2a e ^ t ) 2 — (da — 2a earf ) 2 J . The Lagrangian L sc gives rise to the low 
energy gapless mode - named Goldstone mode. The presence of the Coulomb interaction 
makes the Goldstone mode massive. This can be seen following the method presented in 
ref 12. The Coulomb interaction is replaced by a Hubbard- Stratonovich field 5o, which play 
the role of the temporal gauge field. L sc+C ~ \p ^(d t a — 2a ) 2 — (da) 2 + -^r^\q\ 2 a a 
"e" is the electric charge and |g| is the momentum. Integration of the a field changes the 
Goldstone mode into a massive mode! 

The d x 2_ y 2 case in 2+1 is different! In this case the paramagnetic response is not zero and 
gives rise to a dissipative mode which does'nt become massive in the presence of the Coulomb 
interaction. The difference is caused by the polarization diagram (see ref. 10, 12). Instead 
of rf ' 0) ~ p s and if ,j) = we have a branch-cut, Y^ v) (q) = —}=(-q 2 g^ + q»q v ). As a 

result we have; L d ~ wave ~ \(d^a - 2a^ ! % ul T^' u) (q, uj)(d v a - 2a£ rt )_ f _ w - \p(da - 2a ext ) 2 . 
The most important difference in L d ~ wave [ s the fact that the term (d t a — 2a'^ t ) 2 is absent, 
instead we have the branch-cut terms with the dissipative behavior! In order to consider 
the Coulomb interaction we use again the Hubbard Stratonovich field, a M = <5 M . (Zo- 

T' d ~ wave ~ \(3,a - 2aXA^\Q^)(d,a - 2a„)_ f ,_ w - l -p(3a) 2 + ~a ~a 
z z e z o7r 



Now the situation is different, instead Y^'^ ~ p and Y^ ~ we have; Y^ ^ ~ 2 f ^ 
and ri 8 ' 8 ' ~ — ; = As a result, dimensional analysis shows that the term 7r?\q\ 2 aoan 

2 SyZ-^ + lql 2 ' J 2e2 1 U U 

is negligible with respect the term if ' ^, uj)a a . For this reason we will ignore for the 
remaining part the effect of the Coulomb interaction. 

Next we present our derivation. Our starting point is the nearest neighbor pairing action 
for a two-dimensional square lattice. We start with the superconductor action S and the 
partition function Z. 



Z = J D^D^*DC\DC^DC\DC i exp(iS) 
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H=x,y 



[$(r,r + ^;t)(C|(r,t)C;(r + ^,t)-q t (r,t)q t (r + ^,t)) + / i .c. 



H=x,y 



$(r,r + d fl ;t)\ 2 



(1) 



The action S is characterized by the nearest neighbor pairing fields $ and $* with the 
coupling constant A and the external field a^f*. We will perform a set of transformations 
which leave the partition function invariant. We start by replacing the Fermion field C a (r) 
in terms of chiral fields R Tya and L r ^ ( the right and left chiral fermions). r = 1 corresponds 
to the nodal liquid with Kp = K T= i = anc ^ r = 2 for Kp = K T=2 = f~) 

(see Ref. [5, 6]). K T= \ corresponds to the Fermi surface with the normal in the direction 
<=i = (x + y)/ \/2 and r = 2 is rotated by 90° into the direction e 2 = (x — y)/ \/2. We 
replace the Fermion field by C a (r, t) = J2 r =i 2 [exp(ii? T -r)R T ^(r, t) +exp(— iK T ■f)L T:(7 (r, t)]. 
Following Ref. [11] we have parametrized the superconductor order parameter $ in terms 
of the amplitude p(r) and the "Nambu Goldstone" phase a(r). This leads to a change of 
the measure from DQDQ* to DaDp. We take the saddle point of the action S and find 
that the minimum occurs for p(r,r + d x ) = —p(r,r + d y ) which corresponds to the d x 2_ y 2 
symmetry. We obtain for the saddle point, A = 2y / 2|p| = (f ) 2 f ■ In agreement with Ref. [5] 
we introduce the spinors ^\ = {IV T ^,L\.^) and x\ = (^f i) w ^ n t — 1,2. In terms of 
these spinors we obtain two Dirac Fermion representations with the Cartesian axes rotated. 

Next we perform a gauge transformation [6] which replaces the spinor fields ip T and Xr 
by the neutral nodal fermions ip T and Xt- 




Xr(r,t) = exp(-i—a(r,t))x T (r,t) 



(2a) 



$(r, r + d^; t) = exp(-(ar, t))p(r, r + d^; t) exp[-a(r + d^, t)] 



(2b) 
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From eq. 2a we learn that the spinor field is given as a product of the neutral spinor and 
the boson field exp(e _WJ3 |(r, t)). In terms of the fermions fields we have i? TjCr = e % ^R Tt(7 and 
L T (J = e l ^L T ^ where e ia ^ 2 carries the charge of a half Cooper pair. 

As a result of the gauge transformation given by Eq. (2) the action S is replaced by S. 
We express the action S in terms of the spinors ip T , ip T = V47°> Xr, and Xr = Xr7°> where 
7° is the Dirac gamma matrix. The Dirac action is obtained with the help of the derivative 
expansion. We obtain in the continuum limit two parts: The first part is linear in the gauge 
fields (d^a — 2a^f*) and is given by two Dirac actions Sq T \ t = 1,2 (the two directions see 
refs 5,6). The second part is proportional to (d^a — 2a e * t ) 2 and represents the diamagnetic 
action, S^ a . 

The Dirac action is characterized by two velocities, v i 5 i = v 2)2 = 1 and f lj2 = v 2t i = A = 
(The first index corresponds to the nodal liquid r = 1, 2 and the second one to the 
direction fi = 1, 2). We take the expectation value of (.S^with respect to the free Fermion 
and generate a diamagnetic term [see Eq. (4)]. The diamagnetic mass in Eq. (4) is given by 
r (A) = Af , r = ^ f n 27r f- . cos2g o , where "A" is the ultraviolet cutoff. The explicit 
calculation of the diamagnetic term is given in ref. 6. See in particular eqs. 14, 19, 27 and 
28 in ref. 6. 

For the remaining part we restrict the calculation to A ^ and as a result we obtain two 
Dirac equations coupled to gauge fields in 2+1 dimensions: S = Sq + S T dia , t = 1,2, 

Sq= J d d x[4) T {ffi -m + A T )i) T + x T (iQ - m + A t ) Xt ], m^O (3) 

S T dia = J d d x{- r -^(A T , ) 2 }, d d x = d 2 xdt (4) 

where A T ^ = ^A T ^, d = v T) ^d^ v T ^ = (l,v Ttl ,v Tt2 ), fi,i = v 2 ,2 = 1, v 1>2 = v 2 ,i = A. The 
form of Eqs. (3) and (4) has been obtained after performing a derivative expansion in the 
action in eq. 1. A 2 fi in eq. 4 is the spatial component of the gauge invariant field A^. 

Eq. (3) represents the neutral Dirac equation coupled to a gauge field A T which is a 
gauge invariant quantity. Therefore, the neutral particles ip T and Xr couples to the gauge 
invariant field, Aq. The gauge invariant form follows from the combination of the Nambu 
Goldstone phase and the vector potential a^f*. 
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A-=i,o 

A T =2,0 

A T=1 o 



-d 2 a) 



a 2 



-,ext 



= A T=2 ,2 = a = (a e xt - -d t a) 
= A T=1 o = 



(5) 



It is important to mention that the gauge field in eq. 3 is similar to the gauge fields in 
1+1 dimensions [9] the reason being that the second component is zero A T=12 = A T=2jl = 0, 
(the index r = 1, 2 corresponds to the Fermi surface and the second index represents the 
cartesian direction). Combining the two gauge fields A TtfJi , r = 1,2, we obtain a gauge field 
in 2+1 dimensions, A T= i t o = ai, A T=2> o = a 2 , A T= i t \ = A T=2j2 = a^. The diamagnetic term in 
eq. 4 give rise to a Meisner effect, controlled by the penetration depth, (r (A)) _1 evaluated 
in eqs. 27, 28 in ref. 6. 

Next we will compute the single particle Green's function. We will consider first the 
neutral part. 

G T=1 (r,t) = ((TVv=iM)^ =1 (0,0)>> 



TV 



DA T=1 DA T=2 ^ + 7 ^=i)- t 1 ; o,o e 



°dia 



c(t=2) 



■(det(id + A T=1 )det(ffi + A T=2 )) 2 e lb ^ 



(6) 



The presence of 7 in eq. 6 is due to the fact that we compute %l) T {2)ip\.{l) = ^(2)^(1)7° 
and not ^ T (2)^j T (l). "N" represents the normalization constant. In order to find the single 
particle Green's function in eq. 6 we have to evaluate the determinant; det(i<5 + A T ) = 
expiS eff (A T ). 



W^=lny£- Tr 



n=l 



16 — m 



:-iA T ) 



S 2 (A T ) + Sj(A T ) 



where 



S 2 (A T ) = - I d d qA T ^q)T^(q)A T 4-q) 



with 



-pflU 

1 2 



(q) = 2 lim [ 



d d klr 



7 M (A; + m)Y{k + q + m) 
{k 2 — m 2 )[{k + q) 2 — m 2 ] 



8^q~ 2 



'-q 2 g^ + q»q v 



(7) 



(8) 



(9) 



Si (At) represents the non-Gaussian part. In eq. (9) we used the notation, g^ v = 
for fi ^ v and g 00 = g n = 1. For the momentum q we use g M = (q°,q), q° = OJ and 
q 2 = q^q^ = (q ) 2 — (q) 2 . In eq. (9) we observe that the "branch-cut" for oo/\q\ > 1 gives rise 
to an imaginary polarization [10]. This describes the threshold for destroying the electron 
pairs and creating normal electrons. 

We combine Y 2 w (q) given in eq. 7 with the diamagnetic part S^ a and find the low energy 
action in the gauge invariant form. 



S 2 ,eff(A lj 



d d qA,(q)t^(q)A,(-q) 



(10) 



where 



If %) = if %) + lf 2 >(,), T 2 (q) = t£>%) = t^\q) = if %) - r (A) 
^ 2) (q) = lf%), t^\q)=T^\q) (11) 

In eqs. 10, 11 f ^ ,v \q) is the vertex for the fields a M expressed in terms of the original vertex 
if ^(g) (For the fields A Ttf ,see eqs 8,9). 

Using eq. 10 we obtain the "photon" Green's function. 



D 00 (q,uj) = (a (q,u)a (-q,-u)) 

UJ 2 ~ (Aq 2 ) 2 



= -iA 



+ 



uj 2 - (A Ql y 



-or 



q 2 + (Aq 2 ) 2 -ie) 1 / 2 8(-u 2 + q 2 + (Aq ± ) 2 - ie) 1 / 2 



A-»l 



-tA8(-u 2 + \qf-iefl 2 
2uj 2 - \q] 2 



(12) 



D n (q,u) = (ai(q,u)ai(-q,-u)) 



= iqlA 



k/1 



q\ + (Aq 2 f 



A-t-1 



i8A 



52 



l(-u 2 + q 2 + (Aq 2 ) 2 - ie) 1 / 2 

ql + (A qi ) 2 
8(-u 2 + q 2 + (A gi ) 2 -iey/ 2 

(-u 2 + \q\ 2 -it) 1 ' 2 



ro(A) 
-ro(A) 



\q\J \q\ 2 -8r Q (A)(-u 2 + \q\ 2 -ief' 2 



(13) 



D 22 (q,u) = (a 2 (q,uj)a 2 (-q,-uj)) 



A->1 



tA8[^ 



(-u 2 + \q\ 2 -ief/ 2 



\q\J \qV-8r G (A)(-uj 2 + \qA 2 -ity/ 2 



(14) 
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In the absence of the gauge field A T the neutral fermion Green's function is given by 
G° T=1 {q,uj) = ~?r|~p3p~^ ■ Using the representation in eq. 6 and the dissipative photon 
Green's functions given in eqs 12-14, we compute the fermion self energy, to the one loop 
approximation we find that the self energy is given by Yl T (q,u). 



■<r=l 



r A d 2 k r 

J (2k)*J 



D n (k, n)+D 00 (k, Q) 



(15) 



2vr (Q - uf - (k - q) 2 

The scattering rate in the vicinity of the nodal points, \q\ ~ is given by the imaginary 
part of eq. 15. The main contribution is due to the D 00 (k, Q) photon. Performing a contour 
integral with the pole at Q = u + \ k — q\ ~ uj + \ k\, we find in the limit of A — > 1 the result; 



7T J \ 1 + y 71 A A A 

(16) 

Eq. 16 has been obtained in the limit A — > 1. The reason for this being that according to 
eq. 12 the calculation becomes simple in this limit. Eq. 16 represents the "u;" dependence 
of the scattering rate at T = 0. This result is much stronger than the oo z result known in 
the literature [3]. 

At finite temperatures T < A eq. 16 gives for the scattering rate, ^ ~ ^T^Xx)- Moving 
away from the nodal point \q\ ~ gives a reduced scattering rate for \q\ > uj, 
t ~ ^Kf ) 2 - (x) 2 ]- These results explain the recent ARPES [1] and conductivity [2] 
experiments. 

Next we consider the Green's function for the quasi-particles, 

G T (r,t) = «7V T (r,f)#(0)» ~ (e^WM^t) (17) 

The phase a(r,t) is given in terms of the gauge field a M (for ajf* = 0). 

Integration of eqs 5 with the condition a e * f gives for the Nambu Goldstone phase the 
result: a(x, y, t) — 2[ J* ao(0, 0; t')dt' + J* ai(x', 0; t)dx' + ci2(x, y'; t)dy'\. Using this rep- 
resentation we compute C(r;t) = (exp(— ^(a(r, t) — a(0;0))) a . For simplicity we evaluate 
the equal time correlation C(R = R x , R y ; 0), 
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C(R = R x , R y ; 0) « exp[--((a(x + J R x ,y + i2 I/ ;t) - a(a;, y; t) 2 ) 2 )] 

r 1 /* . fa fn v 1 - cos 9x-Rx n ,1 + cos^iC 
= ex P[-^ / / dg[Dn(g,w) 5 D 22 (q,uj) 

^ J J 1 X 1y 

^ exp[-Ar (A)i? 2 ]; R = J (*Z)2 + Q)2 {18 ) 

"a" is the lattice constant and "r (A)" is the dimensionless inverse penetration length. 
Due to the fact that the correlation C(R) is short range allows us to replace the Green's 
function G T with the neutral one G T (r,t), G T (r,t) ~ G T (r,t). This replacement is justified 
for distances r < ^( ^(A) ) 1 ^ 2 ) sma ^ er than the pair correlation function. 

Next we want to comment on the effect of the nonlinear term Si(A T ) (given in eq. 7). 
Dimensional analysis suggests that our problem in 2+1 dimensions is equivalent to Q.E.D. 
at 3+1 dimensions. This shift in dimensionality follows from eqs. 8,9. To the one loop 
approximation the nonlinear terms Si(A T ) renormalizes [13,14] the vertex T 2 u (q) (see eq. 3) 
to V 2 v {q) ~ |[1 + ^2 h^z^r)] 1 / 2 ' nfX?)- Therefore the use of the normalized vertex f 2 u (q) 
will not change significantly the result in eq. 11. 

In the last part we consider the effect of the Coulomb interaction and show that it can 
be ignored. Following ref. 12 we replace the Coulomb term by a Hubbard field a with the 
action ^2 ^r5 (<f, a;)a ( — q, — cv). a plays the role of the scalar gauge field. We replace in 
eq. 5, a ext = 0, a[f' = a and find for the Nambu Goldstone phase "a" couples to the scalar 
field d the Lagrangian L(a, a ) 



L(a,a ) = a(q,u) uj 2 f ( 2 ' 0) (q, u) + |g1 2 r 2 (<f, u) + ujq^' 1 ' (q, u) + Lvq 2 r^>(q, u) 



(0,1), 



=1(0,2). 



a{-q,-u) 



+ -ao(q,w) 



An 



a (-q, -u) 



+ ^ a (<7» [e(^)r^°' 0) (g,w) - iq^'^i^uj) - iq 2 T ( 2 v ' z> {q, w)J a (-q,-u) +/i.c.(19) 
The explicit form of the parameters in eq. 19 are obtained from eq. 11. In the limit A — > 1 



^(o- 2 )^ 



we obtain: if %~>) = (2cu 2 - \q\ 2 h(q,uj), T 2 (q,u) = \ff-y(q,u>) - r (A); T^>(q,uj) = 
uqiifau), % — 1,2; 7(9,0;) = |(— cj 2 + \q\ 2 )~ l l 2 . From eq. 19 we observe that contrary to 
the "s" wave case (ref 12) the vertex e 2 f(°'°)(<f, to) + ^7^5 — not a constant (see ref. 
12). Therefore the excitations will remain massless in the presence of Coulomb interactions. 
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In order to obtain the Nambu-Goldstone excitations we integrate the field ao and find: 



L(a) = -a(q,u) { \qf 



... , f(°-°)(g;a;)(e 2 47r)- 1 -7 2 (g;a;) 

w Ainin/-. ; + 7(9,^) 



r(°>°)(g» + 



It! 

e 2 47r 



- r (A) + |gl 2 7(g, w)] } -w) 



In the limit \q\ — > we look for massless solutions a; = We substitute a; = into the 
last equation and find a polynomial P(z) for "z"; 

= z 2 (2z 2 - 1)(1 + (e^Tr)" 1 ^ - z 2 )V2) - z 2 - r (A)(l - z 2 ) 1 / 2 ^ 2 - 1) = 0. Massless 
solution exists for z = z(r ) ~ \J ^Jlly r^- [■ We find finite solutions for z ~ ^(^o) ( re& l an d 
imaginary.), therefore the excitations are massless in the presence of Coulomb interactions. 

To conclude, a new explanation for the scattering rate is presented. Instead of the 
phenomenological explanation given in refs 15, 8, we show that the neutral quasi-particles 
are scattered by the dissipative collective excitations giving rise to a linear temperature 
damping at the nodal points, in agreement with the ARPES [1] and conductivity [2,3] 
experiments. The absence of the linear scattering rate in the YBCO superconductor might 
be due to the absence of the dissipative mode caused by the hopping in the "c' direction. 
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